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Abstract. We address the question for existence and uniqueness for the biharmonic equa- 
tion on Lipschitz domains. In particular for the Dirichlet biharmonic problem on 13 C M", 
we show solvability for data in L^, 2— < p < 2{n — l)/(n — 3)+. This result complements 
known counterexamples due to Mazy'a-Nazarov-Plamenevskii and Verchota-Pipher, and is 
thus sharp at least in dimensions four and five. 



1. Introduction 

In this paper we study the problem for solvabihty of the biharmonic equation with 



data on a Lipschitz domain D C 



(1) 



'', n > 4. For the Dirichlet problem 
A^n = on D 



U\dD 

du , 



/ 



\dD 



on 
on 

9 



dD 



M{Vu) e LP{dD) 

we show that there exists an unique solution as long as {f,g) G L\{dD) x U'[dD) and 

2 — e{D) < p < 2{n — — 3) + s{D). This result settles in positive a question posed in 

[17] and should be viewed as a natural extension of the three dimensional results in [15]. 

For the Laplace's equation, which is the standard threshold for elliptic boundary value 
problems, the problem is far better understood. By the results in [6], [22] 



(2) 



Am = on D 

u\q£, = f on dD 
M{u) G LP{dD) 



is uniquely solvable as long as 2—e{D) < p < oo, / G W^dD). Moreover, one has appropriate 
estimates of the solution in terms of the data 

The important endpoint p = oois included, because of the validity of the maximum principle. 
For the regularity problem. 



(3) 



Am = on D 
u\dD = f on dD 
M{Vu) G LP{dD) 



one has solvability and uniqueness provided 1 < p < 2 + e. Moreover, 
(4) I|M(Vm)||^,(,^) < 



2 



STEFANOV, VERCHOTA 



Even though p = 1 cannot be included in the range of (4) one has an appropriate replacement 
if we restrict the data / in the atomic Hardy spaces with one derivative Hl{dD). Then 

(5) l|M(Vtz)||^.(,^)< 11/11^.(9^). 

In a recent work [19] , we have been able to extend (5) in two dimensions to the sharp estimate 

ll^(V«)IL./3(aD)^ll/llj^V3(g^). 

The methods of [6] clearly showed the importance of having well-localized data such as 
Hardy spaces' atoms. Not surprisingly, one proves (4) by interpolating between (5) and 
the estimates of [22]. The results for the Dirichlet solvability can be obtained as a dual 
statements to the Rp results. 

For higher order equations however, it was not immediately clear that solvability would 
follow from solvability for atomic data as it turned out in the harmonic case. In fact, for 
the regularity hiharmonic problem (see Section 3 below) a counterexample due to Mazya- 
Nazarov-Plamenevskii shows that one cannot expect solvabihty for data in U'[dD), p < 4/3 
if C R", n > 5 . This in particular prevents estimates in for the solutions corresponding 
to atomic data. Subsequently, Piphcr and Vcrchota [15] showed that for four dimensional 
domains, one cannot solve the regularity problem uniquely unless p > 6/5 and by the same 
token estimates for the solution corresponding to atomic data necessarily fail. 

In dimension three however, Pipher and Verchota [15], basically carried out the atomic 
approach. They have showed appriori estimates for solution of the biharmonic regular- 
ity problem with atomic data. As a consequence, they were able to obtain existence and 
uniqueness for the Dirichlet problem, together with the estimates 

ll^(VM)|Lp(a^) < WfllLlidD) + \\9\\LP(dD)- 

for 2— < p < oo. In [16], by using the atomic estimates for the Green's function they have 
shown a maximum principle, together with a solvability in the Lipschitz class C^. In [17], 
the above techniques have been further developed to show solvability for higher order 
elliptic operators and an appropriate maximum principle in the three dimensional case. 

Dahlberg and Kenig [5] used similar approach for the related three dimensional Lame 
system to show maximum principle and solvability. 

For the stationary Stokes system, Z. Shen ([18]) have shown maximum principle in the 
three dimensional case and some Sobolev-Besov type estimates (with a derivative loss) in 
the higher dimensional case. 

Evidently, there must be some obstacle to prove estimates in the higher dimensional case 
(n > 4). While there are counterexamples for the biharmonic equation showing such esti- 
mates must fail, one has neither proof nor a counterexample for the solvability of the 
Lame and the Stokes sytems in dimensions higher than three, when p is away from 2. On 
the other hand, one should point out that the theory for all of the problems mentioned 
above has been developed in all dimensions ([17], [4], [7]). 

The purpose of this paper is to shed some light on how to obtain sharp If estimates, 
when p is away from 1. To this end, we still need to exploit the basic idea that whenever the 
data is compactly supported the solution and its derivatives somewhat decay (on average) 
away from the support. The underlying difficulty with this approach is that for general 
data one does not have any sort of reasonable decomposition into atoms. That is why we 
measure the solution in a new familly of weighted U' norms with weights acting on arbitrary 
scales to accomodate various profiles of the initial data. We develop the corresponding (real) 
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interpolation theory for these spaces, so that our LP estimates follow from estimates in LP' 
and H^. 

We believe that appropriate analogues of these spaces can be used to obtain estimates 
for the other elliptic boundary value problems mentioned above. We hope to report on these 
questions in a later paper. 

Acknowledgements The authors would like to thank Jill Pipher for numerous conversations 
on the topic and for pointing out some of the references. 

The organization of the paper is as follows. First, we give some background material for 
harmonic functions. In Section 3 we outline the LP theory for the Dirichlet and regularity 
biharmonic problems. The main estimates are in Section 5 followed by the definition of the 
weighted LP spaces and their real interpolation properties. We state and prove our main 
results in Section 8. Finally, we offer some conjectures for the open problems alluded to 
above. 



2. Some preliminaries 

For simplicity, throughout this paper we will restrict our attention to Lipschitz domains 
D above graphs, i.e. for a fixed Lipschitz function with compact support Lp 

D^{{x,y):y>^{x)}^W. 

It is clear that the general case of non compactly supported ip can be obtained with the 
usual approximation techniques. We will always consider dD as being equipped with the 
surface measure da = -^1 + \S/ip\'^dx. Denote the Lipschitz character of the domain D as 
L — ||V<^||^. We will frequently use the non-tangential boundary cone V{Q) C D associated 
to every point Q e dD. 

r(Q) = {Y e D ■.\Y -Q\<{1 + L/10)dist(r, dD)}. 

The non-tangential maximal function with respect to V{Q) of a function m : D — > C is 

M{u){Q)^ sup \u{x)\. 

xeriQ) 

Define the tangent vectors and tangent derivatives by 

^ F = {Tj, VF) = DjF + ^DnF. 



9Tj ^ ' ^ ■ dxj 



Let ujn be the surface area of ^. Then 

|2-n 



G{x) = -^^^—- (n>2) 

is the fundamental solution for the Laplace's equation in M". Define also the single and 
double layer potentials S and /C by 

S(f)(X) = p.v. J G{X-Q)f(Q)da{Q), xeR"\dD 

dD 

/dC 
-^{X-Q)f{Q)daiQ), xeW^\dD 



dD 
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We also define the formal adjoint of K, 

/f)C 
-^{X-Q)f{Q)da{Q). 

dD 

Prom the LP boundedness of the Cauchy integral on Lipschitz curves [3], we have for 1 < 
p < oo 

(6) ||M(/C/)|| 

LP{dD) ~ ll/llLP(aD)' 

ll^(V5/)llL.(az.)<ll/llLnaD) 

Based on (6), the usual density argument and the jump relations, one defines the singular 
integral operators K^, K_, S^, acting on the boundary dD as 

/C_/(g) = hm ICfiQ) = -IfiQ) + JCfiQ) a.e. 
5+/(g)= hm SfiQ) 

X^Q;XGr{Q) 

S.fiQ) = lim SfiQ) 

The following theorem is essentially a reformulation of the existence and uniqueness state- 
ments for the regularity and Dirichlet problems (see Theorem 2.4 in [15]). 

Theorem 1. (Dahlberg-Kenig, Verchota) There exists e — £(L) > 0, so that 

K:+ : I/idD) ^ I/idD), 2-£<p<oo, 
/CI : L^idD) L'^idD), l<q<2 + e 

are invertible mappings. Moreover, 

(7) ll^(V5/)||^,(,^) ^ \Kf\LiaD) - mLHany 

At this point, one is tempted to say that VSf is an invertible operator, which is heuris- 
tically the case in view of (7). To make this statement precise, recall that one usually uses 
the "preferred" direction X„ to obtain 

\\MiVSf)\\^ - \\MiDr,Sf)\l (Stein's Lemma). 

Theorem 2.7 in [15] gives invertibility of DnS± on Li^idD). In particular DnS± is one-to- 
one map in L'^idD). By (7) and the density of L'^idD) in L^idD), it follows that DnS± is 
one-to-one map in U'idD) as well. Thus, 

Dr,S± : L^idD) U'idD), 1< p < 2 + £ 

is an invertible operator. 

We also remark, that the (small) numbers e that will appear frequently in our discussion, 
will not be the same at every appearance (although one can surely take the smallest one that 
appears and state the theorems with it). That is why, we will sometimes enjoy the liberty 
to denote by A— a number which is equal to A — e for some potentially small £ > 0. 
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3. THEORY 

We consider the Dirichlet and regularity problems separately, partly due to the technical 
issues and ambiguities arising in the definition of the regularity problem. Ideally, the regu- 
larity problem would ask for a biharmonic function u with a prescribed DnU and some second 
derivative on the boundary dD. However, even formally one cannot define two derivatives 
on the boundary due to the smoothness restrictions on the function (f. 

3.1. Dirichlet problem. For the Dirichlet problem (1) we have 

Theorem 2. (Dahlberg-Kenig-Verchota,[4\) There exists e > 0, such that (1) is uniquely 
solvable whenever 2 — e<p<2 + e and (/, g) e {L\{dD) x U'{dD)). Moreover, one has 

(8) ll^(Vti)|Lp < \\^u\aD\\Lv{^BD) ~ {WlWhlidD) + \\g\\LP{dD))^ 

(9) \^u{x)\ < dist(x, 9r>)-("-^)/f . 

Even though Theorem 2 is important in its own right, we would like to somehow relate the 
existence and uniqueness statement in it to the invertibility in U'{dD) of certain singular 
integral operator. We will then essentially follow the approach from Theorem 1, to reduce 
the question for solvability of the regularity problem in to the invertibility of the adjoint 
operator in W {dD). The way we choose to set up the regularity problem will be largely 
dictated by our goal to have the invertibility of the adjoint operator essentially equivalent 
to the solvability for the regularity problem. 

This program has been carried on in the n dimensional case in [15]. Let B be the fun- 
damental solution of the bilaplacian, A^B = S{X). According to Section 2 in [15], one can 
represent the solution as 

G{X-Y)-^)Cf{Y + ten)dY, 

D 

where is the unit vector in the Xn direction. Then the divergence theorem yields 

n „ 

(10) Diu(X) = J {N^Dk-Nl,Dk)DkDiB(X-P))C+f(P)dP + 

+ j {N^Dk - N^Dj)DnDiB{X - P)dP + 

dD 

j j N^gD^GiQ - P)fiQ)dQdP 



+ 



dD an 



for X e R" \ dD. 

Lemma 3.6 and the proof of Theorem 3.7 in [15] show that to prove solvability for the 
Dirichlet problem in i/, it will suffice to show that 

T : / ^ D^u\an 

is invertible, where Dnu\do is the singular integral operator in (10) corresponding to i = n. 
Since the invertibility of T on follows from the Rellich identities (see (3.4) in [15]), 
Theorem 2 for p = 2 follows. 

Remcirk The extension to2 — £<p<2-|-£is automatic due to the Calderon's method 
[2]. In fact, it has been recently been shown in [9] (in a much more general situation) that 
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the set {p\T : U'{dD) L^[dD) is invertible} must be open. In particular, since we have 
verified that 2 is in the set, one has invertibihty on a whole interval 2 — e<p<2 + e. 

3.2. regularity problem. Throughout the paper, we will consider the regularity prob- 
lem 

The following theorem states that this particular version of the regularity problem has unique 
solution, obeying the usual estimates away from the boundary. 

Theorem 3. (Pipher-Verchota, [15]) There exists e = e{L) > 0, so that Rp is uniquely (up 
to a linear function) solvable for {f,g) G {L\{dD) x Wi^dD)) and 2 — e<p<2 + e. In 

addition, there are the estimates 

• |VVm(X)| < dist(X,(9D)-("-i)/p, 

. ||M(VV«)||^,,^) < E,(||Vt,/L,(,^) + \\9h.i9D))- 

4. Solvability of Rp is equivalent to the invertibility of T* 

In this section, we actually prove the equivalence of the Rp solvability 1 < p < 2 and the 
invertibility of T* on U'[dD). We remark that although these results hold true for the full 
range 1 < p < 2, we will really need them only for the range 2{n — l)/(n + 1)— < p < 2. 

4.1. Invertibility of T* implies solvability of Rp. Let h he a harmonic function, such 

that M(V/i) e L'^idD). Note that \Vh{X)\ < dist{X, dDY^-""^/^ and 

\V^h{X)\ < dist(X,(9D)(-i-")/2 as dist{X,dD) oo. Take to > 2max.Rn-np{x) and Xq G 
R"~^. Following [15], one defines a primitive function if of /i by 

t oo 

(11) /f(M)^/M...M.-/(M..^)-M.o. 

to to 

Based on the properties of h, it is not difficult to check that for n > 4, the function H 
is well defined and \V^H{X)\ < dist(X, (9D)(i-'^)/2 and M{V'^H) e L'^{dD). For data 
{9o,gi) G i.L\{dD) X L'^iydD), take h = }C{go) — S{gi). Fix X° = (xo,so) ^ D, where xq was 
chosen before. Let / G L'^{dD). With H defined by (11) set 

D 

= H{X)-G{D,,Sf){X). 

One can check (cf. Section 4, [15]) that u is a solution to R2 as long as one can select 
/ e L'^idD) so that 

(13) T*f = ^D,,G{D^Sf) = ^ICgo - ^S.g, e L\dD) 

Therefore, the invertibility of T* on implies solvability for R2. One obtains similar 
statement for any p : 1 < p < 2, i.e. if T* is invertible on LP{dD), then Rp can be solved 
uniquely. 
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4.2. Solvability of Rp implies invertibility of T*. We show that solvabihty of Rp implies 
the invertibility of the operator T* : U'{dD) — > U'{dD), thus making these two statements 
equivalent. The argument is essentially a reprise of the proof of Lemma 6.2 in [15]. 

Suppose Rp is solvable for some p in the sense of Theorem 3. Consider the "reduced" 
regularity problem with zero DnU data 



(14) 



A^M = on D 

Dnu\dD = on dD 



MiV^u) e LP{dD) 



a on 



dD 



\LP{dD) 



in our estimates. Then T 



where a e L'^{dD) f]LP{dD), but we will use only 

is well defined. Define the harmonic function h — —S{a) and H by (11). Note also that 
according to (12) and (13) the solution can be written as 

Ua^H -G{DnS{T*-^a)). 

In particular it follows that Aua = DnS{T*~^a). Hence, by Theorem 1, Stein's lemma and 
the assumed Rp solvability, one has 



< 



< 



\\M{VS{T*-^a 



< 



< 



\\M{DnS{T* 



\\M{AUa) 



I < 



\LP{dD)- 



Since this inequality holds on a dense set of LP{dD), one has invertibility of T* on LP[dD). 
Observe that, we have not used our full assumption for solvability of Rp, but only the 
solvability of (14). We will show that the regularity problem (14) has unique (up to a linear 
term) solution in the sequel. 

4.3. Construction of the solution for Rp. We have shown that the invertibility of T* 

on LPi^dD) implies solvability for Rp. Wc show now how to obtain the solutions, if we know 
how to solve the "reduced" regularity problem (14). Hcuristically, the solution to the "full" 
regularity problem is performed by solving a harmonic Dirichlet problem and a "reduced" 
regularity problem. 

We consider the "full" regularity problem 

A^u = on D 



(15) 



on 

Dnu\dD = ./ on dD 



g on 



dD 



Miy^u) e L-p^dD). 

First, define h to be the unique harmonic function with Dirichlet data / and 

ll^(v/i)lL.(az,)<ll/ILnai.)- 

Define a primitive H as in (11). We have 
dh 



dh 



and e LP{dD). Consider now the "reduced" regularity problem with data a — 9~'q^ ^ 
U'{dD). Call the solution u . Then u^H + u satisfy (15) with M{V'^u) e U'{dD). 
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5. Main estimates 



Our main results for the Dirichlet and regularity biharmonic problems are the following. 

Theorem 4. Let D <Z he a Lipschitz domain. Then there is e — £{D) > 0, so that the 
biharmonic Dirichlet problem (1) has an unique solution for 2 — e < p < 2(n— l)/(n — 3) + £. 

Moreover 



Theorem 5. Let D d be a Lipschitz domain. Then there is e = e{D) > 0, so that 
the biharmonic regularity problem (15) has an unique (up to a linear term) solution for 
2{n — l)/{n+ 1) — e < p < 2 + e. Moreover 



The uniqueness statements will be proved in Section 8. Let us only remark that the 
methods are standard and can actually be reconstructed from [15]. We note also that the 
restrictions for p come in naturally in the uniqueness results. Albeit not a proof of sharpness 
of our existence results in dimensions higher than five, this of course gives us some indication 
that might be the case. 

As it was pointed out already, the existence statement for the Dirichlet problem would 
follow from the invertibility of T, when considered as an operator acting on LP{dD) for 
2 — e < p < 2{n — l)/(n — 3)+. Similarly, the regularity problem can be solved based on the 
invertibility of T* : LP{dD) LP{dD), for 2(n - l)/(n + 1)- < p < 2+. Moreover, T is a 
bounded operator on L^^dD) from [3]. Thus, it will suffice to show that T* is invertible in 
the range 2(n — l)/(n + 1)— < p < 2+. However that was a consequence of the solvability 
of the "reduced" regularity problem (14) in the same range. Thus, we aim at solving the 
"reduced" regularity problem in the range 2(n — l)/(n + l)— < p <2+. Observe that in the 
three dimensional case, the lower bound for p is 1 — , i.e. one needs to show that the reduced 
regularity problem has solution, when the data is in H^{dD) ( it is quite standard in these 
type of problems to avoid and consider instead H^). Pipher and Verchota have effectively 
used a Cacciopoli type argument to show indeed that such solutions exist and to prove the 
estimates on ||Af (V^m)!!^!^^^,). That was the content of the Main Lemma [15], p. 941. We 
have 

Lemma 1. Let D C M", n > 4 be a Lipschitz domain above graph. Let a be a function 
supported in the unit ball ofdD, a G L^{dD). Then the unique solution u of the "reduced" 
regularity problem with data a satisfies 



+ 11^11 LP{dD)> 

• |Vm(X)| < dist(X,(9D)(^-'^)/p. 



• ||M(V2m) 11^,(3^) < WfllLPidD) + \\9\\LP{dD)> 

• |V^m(X)| < dist(X,c'L>)(^-"Vp. 



(16) 




2 



{{x,cp{x)):\x\r^2j} 



for some positive e — s{D). 



Remark Due to the lack of enough decay in the Green's function associated with low 
dimensions. Lemma 1 seems to be more complicated for D C R^. We will perform an 
additional argument in Section 9 to show that Lemma 1 holds in that case as well. 
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We assume n > 5. Before we go on to the proof, we will need some technical results. One 
has the following Fatou type theorem for biharmonic functions in Lipschitz domains. The 
version below is taken (with small changes) from Theorem 3.9 in [15] . 

Lemma 2. Let D (Z be a bounded Lipschitz domain. Suppose that A'^u — in D and 
M(Vm) G LP{dD) for some 1 < p < oo. Then UjVu have non-tangential limits a.e. on D 
andM{u) e LP{dD). 

The next lemma is a Cacciopoli type inequality for biharmonic functions on Lipschitz 
domains. It appears as Lemma 5.6 in [15]. 

Lemma 3. Let D C be a domain above Lipschitz graph. Let Qi G Q2 D be bounded 
Lipschitz domains and A^-u = in D with M{W'^u) G L'^{dD). Let also 1 < p < 00 and 
d = dist(ili, D \ fl2). Then there is a constant C , depending only on the Lipschitz constant 
and p, so that 



+ d \\u\\Lr'(dDndn2)\\Mi'^ mLP{dD) + 

The next lemma is a somewhat more sophisticated variant of the usual hiding technique. 

Lemma 4. (Hiding lemma) Let {bk} be a sequence of positive numbers, with at most expo- 
nential rate of growth: bk < A2^''. Assume also for some integer I and e > 

bl<B2-''%{bk-i + ... + bk+if' + l). 

Then there exists e' > and a constant C depending on A, B, I, e, N, so that 

bl < CB2-^''. 

The proof of the lemma is elementary, so we omit the details. 
Proof. (Lemma 1) 

Let u be the unique solution to the (14) guaranteed by the regularity theory. Since, 
M(V^m) G LPiydD), we conclude by Lemma 2 that m, Vm, V^m have non-tangential limits 
and the maximal functions taken over some cones with finite height are in U'[dD). Define 

(17) u{x)^, — \ — / ^^^y^^dy. 

M"-i 

Since u is a convolution of a{y, (p{y)) with the Green's function, we get A^u — a. Also, by 
differentiating the integral above, one obtains the estimates 

1*^1 11*^112' 
1*^1 11*^112' 
1*^1 11^112' 

for large x. Note that since Dnu\dD — 0, one has 

(21) Dju{x, (p{x)) = Q^ix, ^{x)). 



(18) 


\u{x)\ 


< 


(19) 


\Vtu{x)\ 


< 


(20) 


\Vtu{x)\ 


< 
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By Lemma 2, one justifies the following calculation. Take a test function ip e C^(M'^~^) and 
perform two integration by parts to get 

/'^ ^ f dtp 
Aip{x)u{x, (fi{x))dx — — / — — Dju{x, (p{x))dx — 

n—1 „ / n—1 \ 

Yl J ^^^) I 5^(Vr,> VD,«)|9Z5 j dx = {il^.a), 



Sn-l 



where in the second to the last identity above, we have used (21). 

Hence, u{x, ^{x)) is a weak solution to the equation A^v^x) = a{x, (p{x)). But so is u{x). 
Thus, by the Weyl's theorem for uniqueness of harmonic functions, (20) and M(V^ii) G 
L'^{dD), we deduce that u{x,Lp{x) and u{x) differ by at most a linear term. By subtracting 
the linear term, we can assume that in fact u{x, ^p{x)) = u{x). 

We now turn to the proof of (16). Note that the case of j < 3 follows from the regularity 
theory. Assume j > 3. We then dispose of the supremum in the definition of M(V^xi) taken 
over the points inside D that are far away from dD. More specifically, take Tq to be a cone 
pointing upward with vertex at the origin and with a large slope, say 100||(/9'||^. If Q G dD 
then the intersection rof]r{Q) consists of points X with dist{X,dD) > \Q\. Define 

Mi{V\){Q) = sup |V^m(X)| 
xeronr(Q) 

M2{V''u){Q) = sup 

xer(Q)\ro 

We will show that 

(22) Mi(VM(g) < |gr'-^"-')/('-^)||a||2 for large Q, 

which implies (16), when M is replaced with M^. To show the pointwise estimate (22), we 
have by the Dirichlet theory, (17) and fractional integration 

and therefore 

(23) |V^i(X)| < ||a||2dist(X,aL>)-("-^)/('-"). 
Consequently, for all X eTq f]T{Q), we get by (23) and interior estimates 

< ||a||2dist(X,aD)-^-("-i)/M, 

and hence since dist{X,dD) \Q\ 
It remains to show the bound 

(24) J M2{V\f < \\a\\lR^-'-'\ 

{{x,<p{^)):\x\^R} 

for i? > 10 and some positive £ > 0. 

For 1 < r < 2, define the Carleson region above Zr — {{x, (p{x)) : \x\ ~ R} as 

= = {(x, t) : R/r < \x\ < Rr, ip{x) < t < 100Ti?||^'||^} . 
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From the regularity result on Q^, we have 

M2(VVm)2< / |VVMp + J] / |Vr,V«|l 

By (20), we get for Q = {x, (p{x)) e dflr n dD 



since \x\ ~ i?. Hence 



n-l „ II ||2 

a L 



II, A. n 

E y ivT,vz.r< 



Averaging (25) in r e (1,2) (see p. 944, [15]) yields 



(26) y M2(VV«)^<i?-iy|VV«p + -^. 

By Lemma 3 with p = 2 — £, we can choose a domain Q3, such that C Q3 C D, 
dist(Q2, -D \ Jla) ~ -R and 

(27) y"|V^«r < l|V«IL.+(az.nan3)ll^(VM|L.-(az.) + 



^2 

+ ^'i«llL2+(aDn9n3)II^C^'«)|lL2-(9D) + 

+ ^"iVM||^2(f,3)||V'K||^,(^^)+i?-lii||^2(n3)||V'M||^2(n3)- 

By the I? regularity theory 

(28) \\M{y^u)\\^,_^^^^ < ||a||^2-(a^) < M^.^sDy 

The boundary terms are estimated via (18), (19) 

(30) lkllL2+(aDnan3) ~ ^(!!_l+)/2- 

The fundamental theorem of calculus and (18) yield 

Putting together (27), (28), (29), (30), (31) yields 

,2 / Il«ll2 I hU^'^iyins) 



(32) / 



^2 



^^""1 ~ Rin-3+)/2 ^ i?(n-2)/2 ^ 



+ i?-lVl.||^2(n3)||V2i.||^,(^^). 
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For some terms, the trivial estimate 

(33) P"4LHns) ^ R'^'lM'^'^A^an) ^ R'^'^h 
will do. For others, we also have 

/ iv^«r <rIJ M,{v\r + 1 M,{v'uA , 

since it is possible to have ^3 f](ro fl T{Q)) 7^ for some Q e dD. However, by (22) we 
easily bound the contribution from Mi(V^'u) 

(34) J M^{W'uy < R-'-^\\a\\l. 
Combining (26), (32), (33) and (34) gives 

(35) / M,(V\)^ < ^ + i?-^/lV«||,.(^3)||M2(VML.(an3nai.)- 

dflindD 

To estimate Jq^_^^qd M2(V^m)^ based on (35), we shall need to apply the Hiding lemma and 
apply Sobolev embedding. First, choose p — 2 — e and denote its conjugate exponent by 
p' = 2+. Note that 

1 n n — 2 3 

p^p' 2~ ~ 2~ ' 
This is easily checked by setting f{p) = 1/p + n/p' — (n — 2)/2 and verifying that this is a 
monotonically increasing function in (1,2) with /(2) = 3/2. By Holder, Dirichlet theory 
and the Sobolev embedding ^^2,2(1^") M^2n/(n-2),i(lR") we obtain 

J\'^u\'<lj\Vu\A R-/^' I J I v«r' I < 



(n-2)/2n 



< i?lWp'-(n-2)/2||^||^||^2^||^^^^^^ < i?-||a||,||M2(V^u)||,,(,,^^,^^ 

Combining this last estimate with (35) yields 

(36) J M,{y\f < i^ + i?-^/'-||a||f ||M2(V2-^ii'/' 

For i? = 2'^, denote 

Ok — TT-r, 

It is clear that 



I«ll2 



2'll^2(V2tx)||^.(an3naD) 



10 

< ^ hk+i- 



' "2 l=-W 



ON THE BIHARMONIC DIRICHLET PROBLEM 



13 



Thus we rewrite (36) as 



10 



-10 



An application of the Hiding Lemma (Lemma 4) to the sequence {bk} yields 

or equivalently (24) for R = 2^^, whence the general case for non dyadic R follows immediately. 

□ 



6. Some interpolation tools 

This section will provide some background material on interpolation spaces. Although the 
facts are quite standard, the definitions of the spaces that we are about to use differ from 
one source to another. Since the interpolation formulas in the endpoint cases can not be 
easily put into an unified framework, many authors prefered to leave them behind. That is 
why, we felt we needed to present the basic theorems, with an emphasis on the ones which 
we will be using in the sequel. The reader might find it convenient to skip this section at 
first and use it only as a reference later, when need arise. 

We start with the definition of the Tricbel-Lizorkin spaces for which wc follow the expo- 
sition in [21]. Let Vo G S^{R'^) with supp C (-1, 1) and ipoix) = 1 for -1/2 <x < 1/2. 
Take ^{x) = ^o{2x) - ^o{x). CaU ^j{x) = ^{2-^x) for j > 1. Then 

oo 
j=Q 

This defines the Littlewood-Paley projection operators 

Sjf{x) == j f{x - y)ipj{\y\)dy, 

which essentially restrict the Fourier transform of / to the annulus : \^\ ~ 2^}. The 
Triebel-Lizorkin spaces F^^^ are defined as as the set of all functions / with 




For a background material on F^^ and their relation with Lebesgue and Besov spaces one 

may consult [21], p. 169. 
Introduce the spaces F*' 

^p,q,(r) — |/ ■ ||ll'^.7/lli<!|lLP,T-(Kn) < ) 

which will used to describe {Fpl,q^-i Fpl,q2){e,p)- We remind that L^'"' are the Lorentz spaces 
and W'^ = LP, while is the weak U space. The following lemma is due to Triebel (cf. 
Theorem 1, p. 184, [21]). 
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Lemma 5. Let 1 < pi < oo, 1 < p2, q < oo, pi P2- Let 1/p = (1 — 0)/pi + 9/p2- Then 

Remark The proof provided by Triebel in [21] does not explicitly state the case pi — 1, 
which we will need. Actually, one can easily follow the argument in [1] for Besov spaces ( see 
p. 153, (2)), where the case pi = 1 is covered. Alternatively, close inspection of the proof in 
[21] shows that the argument goes through in the case pi = 1, p2 7^ 1, if one uses the formula 

It is also a standard fact (see [8]) that the Hardy spaces can be put in the framework of 
Triebel-Lizorkin spaces, namely — F^2- By the Littlewood-Paley theorem = F°2- By 
Lemma 5 one gets for 1 < r < oo 

(37) {H\Lye,r)^Fj,,2,ir) where l/p ^ {1 - 9) + O/q. 



7. Interpolation space 

In this section, we introduce a familly of function spaces X^(]R") in which we will measure 
our solutions. Heuristically, a is a weight index, while p stands for LP integrability as in 
the usual Lebesgue spaces. We have several objectives. First, we would like an inherent 
connection with LP spaces, i.e. we wish LP to be somehow embedded into this familly of 
spaces. In fact, we will show that Xq(M") ^ LP. Second, we would like this familly to 
have a good "scale" properties. While we cannot quite satisfy that with our construction, 
we still have an almost precise formula for interpolation of a pair of spaces with the real 
interpolation method. Lastly, we would like to be able to translate the estimates in Lemma 
1 into estimates for the solutions in spaces. 

To this end, fix a sequence {flm} C M"^, dense in the unit ball of MP and \am\ < 2. 

For the interpolation scheme, we will represent X^ as a factor space of a sequence space 
and later we will use standard interpolation results to find the intermediate spaces. Let 
be the sequence space consisting of functions (fi^m): fi,m : K" — > C, such that 



II (/z,m) Hyp (Mn) — 

;<0,m>0 



j>0 



< OO. 

LP(R") 



Here, we restrict our attention to 1 < p < 2 and a E J, where J is a finite interval (we will 
use below J = {—{n — 3)/2+, 1+)), although there is no principal difficulty to extend the 
interpolation theory of Y^ spaces beyond that. Define 



C 



l{fi,m)e fl Yr. Yl hm = o] 

I l<P<2,(TeJ l<0,m>0 ) 



Note that Y^ ^ 0, since at least the functions in the Schwartz class S have the 

i<p<2,o-eJ 
required decay properties. 



Claim 1. C is a closed subset in each Y^ ,1 < p < 2,a E J. 
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Proof. Take a convergent in sequence (//^) in C. We need to show that the hmit 
{fi,m) G adds up to zero. Fix j > I. We have 



ifl,2 fl,m)X\x-am\<2i 

l<0,m>0 



as n — > oo. 



LP 



Since |a^| < 2 and ^//^ = 0, we conclude 

^//,m(^)=0 for |x|<2V2. 
This imphes Claim 1, since j was arbitrary. 



□ 



Define Xti 



I") = yj'(]R")/C or equivalently as the space of all functions / : M" — >■ C so that 

/ ^ ] fl,mi {fl,m) £ 



/<0,m>0 



inf 



E 

KO,m>0 



J]/,,^(x)V,(2-'|x-a^|)2^'^ 

j>0 



< oo. 



LP(M") 



Since ■0,- is a partition of unity, we observe that for u = 0, we have 



/ Jl,m . J 



> 



LP- 



i,l,m 



Thus we obtain 

(38) Xl{W) ^ LP(Mr). 

Turning to the interpolation issues, we would really like for the complex interpolation formula 

[X/C,Y/C]e^[X,Y]0/C 

to hold. Although heuristically right, the formula might fail since C does not seem to be a 
complemented subspace in Y^ (see the discussion on p. 120 in [21]) and as far as we know it 
is an open question whether the complementability condition is really necessary. However, 
there is the following real interpolation version, which does not require complementability 
of C. It is due to Petunin ([14], see also [21], p. 120). 

Lemma 6. (Petunin) Let {Xq,Xi) be an interpolation couple, and C G XqH Xi is a closed 
subspace of both Xq and Xi. Then for < 6 < 1, 1 < p < oo 

{Xo/C,Xi/C)(0^p) = (Xo,Xi)(e,p)/C. 

We are now ready to state our main interpolation result. 

Lemma 7. Let 1 < pi < 2, 1 < pi < 2 and < 9 < 1. For any (7i,(T2 G J , define 
(7 = (1 - e)ai + ^(72 e J and 1/p = (1 - e)/pi + e/p2. Then 



(39) 



In other words, the indices 9,p in X^. interpolate like the usual weights and integrability 
indices for weighted LP spaces. 
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Proof. By Petunin's Lemma since C is closed in both Y^^ and Y^^, it will suffice to show 
that 

(40) (yi'/,Fi;^)(,,i)^i7 

under the same restrictions on pi,p2 and 9. Write 

<,^(^) = E^^(2"1^-«-l)2^'^'- 
j 

Then 

Ya' = C/^^«/^)^^)), 

where we use the notation Ij for the sequence space indexed by i. Standard interpolation 
results stipulate that the real interpolation method applied to spaces in the form l^{Aj) 
results in spaces in the form 1^{A), where A — [Ai, >l2)(e,i)- In our situation, by the theorem 
on p. 121 in [21] 

= (4,,(L^n</^)rf:^))),C(^''«/(^)^^)))(^,i) = 

= C(^'n<M^)^^)), L^^«,(x)cix))(,,i)). 

At this point, we would have prefercd to have a complex interpolation method applied to our 
spaces, since weighted spaces do not behave very well under real interpolation. Actually, 
there exist general formulas that describe the (real interpolation) intermediate spaces of 
weighted U' spaces, but they are too complicated. Instead, we will pass to the complex 
interpolation functor by the well known relation (see for example p. 102 in [1]) 

(X,Y\o,i)-^[X,Y]o. 

Thus, since 
we obtain 

iYl[\Yl^\e,i)-^Y^^, 

which establishes our interpolation step. □ 

8. Existence and uniqueness 

We start this section with the existence statement in Theorem 4 and Theorem 5. 

Proof. (Existence) As it was pointed out allready (see the discussion after Theorem 5), it 
suffices to show existence for the "reduced" regularity problem in the range 2{n — l)/{n + 
1)— < p < 2+. As always things are reduced to showing the estimate 

||M(V^ <C\\a\\ 

for smooth data a and a constant C which is independent of a. Set the sublinear operator 
Ta{X) — M{V'^Ua^^■^{\.\)){X) for X e dD, where V'o is the fixed smooth cut-off of (—1, 1). In 
other words, we take the solution that corresponds to a'(/'o(| ■ I) instead of Ua- This is done 
in order to localize the problem to data supported in the unit ball. We will show that 

(41) ll'^'^ll^-{„-3)/2+('5''f^) - 

(42) l|ra||^.^(3^) < C||a||^.(g^), 

with the usual identification of functions on dD with functions on R"~^. 
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Thus by the definition of T, it will suffice to prove (42) for functions a with support in 
the unit ball. To this end, take = and I — 0, /o,o = M^V'^u^^a) in the definition of 
||M(V^iia)||x2^(aD) ^-nd observe that (42) follows from the statement of Lemma 1. 

For (41), let a — — (n — 3)/2 + s/4, where £ > is the positive number guaranteed by 
Lemma 1. We expand a — ^ )^mbm in sums of atoms. Similar argument as the one above 

m 

shows that it will suffice to consider b„iipQ instead of b^. We estimate first the contribution 
of atoms with diameter {suppbm) > 1/10. Set 

9 — ^mbm- 

m: size {suppbm)> 1/10 

Clearly g e f] L°°. By Cauchy-Schwartz and Lemma 1 

j>0 j>0 

For the "small" support atoms, set supp b^ — -6(2;^, r^^). Let Im '■ fm ^ 2'"" < 1/10. Choose 
q — q{rn) such that \aq — Zm\ < '"m/lO- Thus, by the triangle inequality, we have 

" m j>0 

Thus, it will be enough to prove for a fixed "small" atom 

(43) E2^H^(^'«^o.J|LH|.-s(.)h2<-^) ^ 1- 

We have by Cauchy-Schwartz and Lemma 1 (with the appropriate scaling) 

E2^H^(^%o6j|L.(i_^^^,i......) < 

j 

< ^2^('^+(n-3)/2-e/2)2U(n-l)/2||^^^^||^^ < 



which establishes (43) and thus (41). 
We show now that (41) and (42) imply 

(44) \\Ta\\^, < ||a||^„ 

for the range 2(n — l)/(n + 1)— < p < 2+. 

Estimate (44) follows from the next lemma with ai — —{n — 3)/2-|-, (T2 = 1+ and 9 — 
(n — 3)/(n — 1)— for p ~ 2(n — l)/(n + 1)— and then by complex interpolation with the 
theory. 

Lemma 8. Let ai < Q and (T2 > 0, such that (1 — 9)ai + 9a2 — and 1/p = (1 — ^) + 9/2 
and assume that for a sublinear operator T 
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Then 

T : LP+ ^ 

Proof. By the real interpolation method, we get 

T : (i7\L^)(e_i) {Xl^,Xl^)^e,i)- 
According to Lemma 7, one has 

Therefore, 

\\Ta\\L, < \\Ta\\^xi^^^,^Xe,r, ^ ll«ll(i/i,L2)(,,,)- 

By (37) 

(45) M^,< llall^o 

Prom the estimate for T, we have in particular 

(46) \\Ta\\^, < \\a\\,, = ||a||^o^^^^ 

The interpolation of F^^ (i) spaces is in fact very similar to the interpolation for the usual 
Triebel-Lizorkin spaces. Set S > and let ps : 1/ps = {I — 5)/p + 5/2. It is clear that 
p < ps < p + 0{6). By an argument similar to those in [21], p. 185, claim (c), with the 
appropriate replacement of with either L*'^ or L^'^, one gets 

where the last identity is the Littlewood-Paley theorem for L^^ . Also 

Interpolation between (45) and (46) with {S,ps) yields 

□ 

It remains to observe that for data a whose support is inside {x : \x\ < 1/2}, we have aipo = a 
and therefore (44) reads 

||^(V'«a)L.<||a|L., 
for 2{n — 1) / {n + 1)— < p < 2+. But in this last estimate, one can rescale to prove 

72, 



||^(V'^^a)IL.^II«llL., 
for data a having compact support. The usual approximation techniques finish the proof. □ 

For the uniqueness of the regularity problem, we refer the reader to Lemma 6.9 in [15]. 
Although the statements include only the case D (Z M.^, one can check that the higher 
dimensional case follows as well. In fact, the proof is a lot easier, since we are not anymore 
in the endpoint case L^{dD), where the predual space is unavailable. 

For the uniqueness of the Dirichlet problem, we follow Theorem 7.1 in [15]. Suppose 
A^u ^ u\dD ^ 0, du/dN = with M{Vu) e L2("-i)/("-3)+. Denote by D the original 
domain, which is translated by one unit up. Define the domain 

D^{X : X/\X\^ e D} 
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and 

u= |xi^-"M(x/ixn. 

The function u is biharmonic witli zero Dirichlet data. To show m = 0, it suffices by the 
uniqueness results of [4] to show that MiVu) G LF'^{dD). Fix a cone V with vertex at 
(0, 1 + ip{Q)) e dD. Define 



X 



G r 



FoUowing the estimates in [15], we find that 

sup \Vu{X)\{Q) e LP{dD), 

xer{Q)nr 



for all p < oo. 



For the iterated Hardy-Littlewood maximal function Ai2{f) = M{M{f)), there are the 
estimates ([15]) 



sup \VuiX)\{Q) 

xer{Q)\r 



< I {M2ivu){p)r+\p\-'-dP < 

L^+{dD) dD,\P\>l 

(„_3)/(n_l) / \ 2/(n-l) 



< 



|p|n-l+ 



dP 



< OO. 



\P|>1 



Thus, ||M(V'u)||^2+(g£,) < oo and by the uniqueness result of [4], u — 0. The uniqueness part 
of Theorem 5 follows. Note that the proof that we have presented breaks down for exponents 
p < 2(n — l)/(n — 3) — . That seems to indicate that 2(n — l)/(n — 3)— is the sharp exponent 
in dimensions n > 6. 



9. The four dimensional case 

In the four dimensional case, one does not obtain Lemma 1 (at least not directly), due to 
the lack of enough decay of the Green's function. Wc need an additional argument. 

In the proof of Lemma 1, one obtains an estimate for the maximal function away from the 
boundary 

which implies (16) for Mi. For M2 an identical argument as in the proof of Lemma 1 with 
n — A, yields 



(47) 



{{x,ip{x)):\x\^2^} 

Thus we have (47) for M(V^Ma) as well. Hence, following the proof of estimates (41) and 
(42), one establishes 



T: ^X^ 



3/4+) 



T : 



^3/4+- 



where Ta = M (V^ Ua^^). According to Lemma 8 that implies 

T:L^/^-{dD)^L^/^-{dD). 
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By rescaling, we get for the "reduced" regularity problem 



(48) 



||M(VV 



LP{dD) ~ II'^IIlpCSI?) 



4/3- <p<2. 



This enables us to go back to the proof of Lemma 1 and improve on our estimate (47). 
Indeed, in the derivation of (26), we have used the Cacciopoli estimates from Lemma 3 
with p = 2—. However, we have at our disposition (48), so we choose p — 4/3—. A quick 
inspection of the proof shows (just as before), that one gets 



(49) 



J M,(VM^<2(-^-)ia||i.(,^), 



which is Lemma 1 for n — 4. Prom there on, the proof proceeds as in the higher dimensional 
case n > 5. 



10. Some open problems 

In this section, we list some open problems for boundary value equations (or systems) 
that are somewhat related to the biharmonic equation. As we have mentioned already in 
the introduction, the common feature in all of those are the maximum principles and Holder 
solvability in three dimensions, while for dimensions higher than three, only the theory 
has been developed. 



10.1. The Lame system. Let D c 

consider the Dirichlet problem 



be a Lipschitz domain and u = {ui, . . . , «„). We 
Am + Vdiv?/ = 0, 

(Dp) u\9D = f 

M{u) e LP{dD) 



and the traction problem. 



Am + VdivM = 0, 

(Tp) {yu + Vv')N\aD = g 
MiWv) e LP{dD) 

We remark that the traction problem seems to be the right substitute for the regularity 
problem for the Lame system. 

For n — 3, Dahlberg-Kenig [5] have proved (weak) maximum principle in full analogy 
with the harmonic case, by reducing to the traction problem. For the traction problem, they 
have succcsfully used the atomic estimates method of [15], together with the appropriate 
Cacciopoli type inequalities. 

For n > 4, there is no counterexample to a maximum principle, but the best solvability 
result might be for 2— < p < 2{n — l)/(n — 3) + , just as in the biharmonic case. Let us only 
remark that as is well-known the solution u must be a biharmonic vector. 



10.2. (Stationciry) Stokes system. 



Au^ Vp 
divu = 

— * 
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In dimension three, Z. Shcn proved (weak) maximum principle for the Stokes system and 
some Sobolev-Besov type regularity results (with a derivative loss) in n > 4. A natural 
question that arises is about the LP solvability for dimensions higher than three. 

10.3. The polyharmonic equation. We will be highly schematic for the definiton of the 
polyharmonic equation. The reader is refered to [17] for an extensive treatment of these 
higher order boundary value problems. For suitable differential operators on the boundary 
Po, ■ ■ ■ , Pm-i {Pi is of order i) 

A'"m = 

Pou\dD = fo, . . . , Pm-lu\dD — fm-1, 

In [17], Pipher-Verchota have shown that the usual L?' theory for all n, m holds and a (weak) 
maximum principle for n = 3, just as in the harmonic case. The results again made use 
of the atomic estimates that we have alluded to earlier, and we ask whether an analog of 
X'P spaces might be helpful to study the LP solvabihty in dimensions higher than three. In 
particular, a potentially sharp estimate in the form 

might hold for all 2— < p < 2(n — l)/(n — 3)+. 
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